Confinement of J\f = 1 Super Yang-Mills from Supergravity 
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We calculate circular Wilson loop of pure Af = 1 super Yang-Mills from Klebanov-Strassker- 
Tseytlin solution of supergravity and proposed gauge/gravity duality. The calculation is performed 
numerically via searching minimal surface of string worldsheet. It is shown that Wilson loop exhibits 
area law for large radius such that Af = 1 super Yang-Mills is confined at large distance. Meanwhile, 
Wilson loop exhibits logarithmic behavior for small radius and indicates asymptotical freedom of 
Af = 1 super Yang-Mills. 



INTRODUCTION 



(N 

J> . Physicists now believe non-abelian gauge theory exhibits some dramatic behaviors at long distance, such as con- 
7-H ' fmement and mass gap, etc.. It is difficult to study these properties using traditional quantum field theory (QFT) 
method due to its nonpertubative characters. However, the remarkable success of AdS/CFT correspondence reveals 
a practical approach to study strong coupling gauge theory 0, 0, 0- The similar correspondence is believed to exist 
in broad class of theories even with less supersymmetry and non-conformal theories (so-called gauge/gravity duality). 
According to continuously breaking supersymmetry, we can expect to study strong-coupling properties of non-abelian 
gauge theory via dual weak-coupling gravity. So far, several nonsingular supergravity (SUGRA) solution have been 
constructed. They preserve 1/2 or 1/4 of the maximal supersymmetry and are conjectured dual to d = 4, Af = 20 
or Af = ljEISGl su P er Yang-Mills (SYM) theories. In particular, the case Af = 1 super SYM is very interesting, 
■ , because it possesses the common features of non-abelian gauge theory at long distance but is hard to be study using 
usual QFT method. For this SYM, two dual SUGRA solution were found: 

'-^ • The Maldacena-Nuhez (MN) solutionQ: It corresponds to a large (N) number of D5-branes wrapped on a 

J> , supersymmetric two-cycle inside a Calabi-Yau threefold. Meanwhile, 6d SYM theory onto D5-branes is reduced 

to pure d = 4, Af = 1 SU(N) SYM in the IR. 

• The Klebanov-Strassler-Tseytlin (KST) solution [g,Lj: ^ describes the geometry of the warped deformed conifold 
when one places M D3-branes and N fractional D3-branes at the apex of the conifold. It is dual to a certain 
Af = 1 supersymmetric SU(N + M) x SU(M) gauge theory. If M is a multiple of N, then this theory flows to 
SU(N) in the IR, via a chain of duality cascade which reduces the rank of the gauge group by N units at each 
cascade jump. Thus at the end of the duality cascade the gauge theory is effectively pure Af = 1 SU(N) SYM. 

The purpose of this paper is to study confinement of Af = 1 SYM from KST solution. We will give a brief comments 
on MN solution at the end of the paper. 

The natural criterion for confinement of pure Yang-Mills theory is expectation value of large Wilson loop satisfying 
area law@. The evaluation on Wilson loop of Af = 4 SYM from supergravity in AdS space has been studied in 
various aspects. It was proposed m|^[lCj that Wilson loops of the CFT can be described in AdS by 

< W{C) > ~ lim e ~( s *- m \ (1) 

where S is the proper area of a fundamental string world-sheet which lies on the loop C on the boundary of AdS, I 
is the total length of the Wilson loop and $ is the mass of the W boson. The presence of the term ^<i> is to subtract 
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divergence from infinity mass of the W boson. It corresponds to UV divergence in QFT. The evaluation for the 
rectangular Wilson loop was performed in refs. Q, ^| and for simpler circular loop was performed in refs. 0, • 
There are no any surprise results due to conformal invariance. The more non-trivial results is from calculation of 
operator product expansion for the Wilson loop when probed from a distance much larger than the size of the loop[I3, 
and from calculation on Wilson loop correlator 13] . The extension of Eq. (JTJ to KST background is in principle direct. 
However, since KST background is much more complicated than AdS space, we can not expect to obtain any analytic 
results. So that although there is a simple scaling argument that large Wilson loops of JV = 1 SYM resulted by KST 
solution satisfies the area raw|l4j. the direct verification is still lacked so far. In this paper we use numerical tools to 
search minimal surface of string world-sheet in KST background which lies on the Wilson loop on the boundary of 
KST space. In order to simplify the calculation, we consider circular Wilson loop of JV = 1 SYM in Euclidean space. 
The proper area of string world-sheet can be expressed as a function of radius of circular loop. We will evaluate both 
of cases for small and large radius. It can be expected that behavior of Wilson loop exhibits the following conclusions, 

• The confinement of N = 1 SYM at large distance. 

• The asymptotical freedom of JV = 1 SYM at short distance. 

• Phase transition from short distance to large distance. 

This study can also be treated as an examination for duality between SUGRA KST solution and JV = 1 SYM. 

The paper is organized as follows. In Sec|H] we extend to calculation on Wilson loop in AdS space to KST 
background and derive Euler-Lagrangian equation which governs minimal surface of string world-sheet. In Scc lIIII 
We show some numerical results such as curve of Wilson loop vs. radius of the loop. The physical results are obtained 
via those numerical results. A brief summary will be devoted in Sec II VI and several details on numerical evaluation 
are presented in the Appendix. 



II. WILSON LOOP OF TV = 1 SYM FROM KST SOLUTION 



The Wilson loop operator in JV = 1 SYM is 



W(C) = ^TrPe^c A i 



(2) 



where C denotes a closed loop in spacetime, and the trace is taken over the fundamental representation of the gauge 
group. The expectation of the Wilson loop can be calculated by supergravity in terms of Eq. Q directly, with action 
S is given by 



S 



2ira' 



7 / d aV^g, 



(3) 



where a 1 and er 2 parameterize fundamental string (Euclidean) world-sheet, 



g a0 = G MN d a X M d X N , (4) 

is induced metric on the world-sheet. In this paper, we focus our attention on nonsingular KST background such that 
Gmn is given by 



ds\ a = h~ 1/2 (T)dx n dx n + h 1/2 (T)dsl, 
h{r) = a J 
= o/(r) 



, x coth x — 1 . . , ,„ . „ ,,/, 

dx : — 5 ( sinh ( 2x ) ~ 2x) 1/3 

sinh x 



(5) 



where r is the radial parameter of the transverse space, a = 2 2 ^ 3 (g s Na') 2 e 8 / 3 . ds\ is the metric of the deformed 
conifoldlljj, 



dsi 



1 



e 4 / 3 K(r) 



J^Adr 2 + (g 5 ) 2 ) + cosh 2 (I) ^(^) 2 + sinh 2 (I) j^tf? 



where e parameterizes length of the deformed conifold and possesses mass dimension —3/2, K(t) is defined by 

(sinh(2r)-2T) 1 /3 



K(r) 



2 1 /3 s inh- 



(6) 



(7) 
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and the 1-form basis g % parameterize transverse five-dimension space. To simplify complicated calculation we do the 
following notions: 

• Set angle parameters of transverse space as constants such that g l = 0. The Wilson loop in TV = 4 SYM with 
non-constant angle parameters has been considered in ref . 0] . When we extend this consideration to M = 1 
SYM calculation becomes very difficult. 

• To define a new dimensionless radial parameter r — e _T / 3 , < r < f . The boundary of KST space is r = 0. 

• We consider that TV = 1 SYM is defined in Euclidean space, and circular Wilson loop is located at x 1 — x 2 
plane. Since we are aiming at finding a minimal surface, the string world-sheet can be considered as a rotational 
surface. Then we define 

x 1 = y|ae 2/3 /W cos 9, 

x 2 = y|L 2 / 3 /(r)sin0, ( 8 ) 

with a dimensionless function /(r). Consequently the string world-sheet is parameterized by r and 9, 

ds 2 = ^e 4/3 /- 1/2 (r)|[/' 2 (r) + ±I(r)K- 2 (r)]dr 2 + f{r)d9 2 } , (9) 

where prime denotes derivative on r. 
According to the above notions, the action defined in Eq. ((3J) is rewritten to 

S = 4i^ 9sN J dr ' J ~ 1/2 /V/' 2 + G(r), (10) 

with 

1 



G(r) = -I{r)K-'(r). (11) 



r 



From action l|10Jl we obtain the Euler-Lagrangian equation which govern minimal surface of string world-sheet, 

f" \{\nIG)'f - ^/' 3 - j(f' 2 + G) = 0. (12) 

The minimal surface implies the initial conditions of Eq. I|12|) should be 

f(0) = R, /'(0) = 0. (13) 

In order to show that this initial condition is taken care of properly, let us find asymptotical solution of Eq. i|12|) for 
r — > 0. Noticing the expressions, 



=> G(r) ^--(1 + 12 In r), (14) 



J W ™ -1^173^(1 + 12 lnr) + 0(r 10 ) 

K(r) ^ 2 1 / 3 r + (9(r 7 ) 

at limit r — > 0, the Euler-Lagrangian equation 1)12(1 approaches to 

//" - (~ + -j-)ff + I Inr + 1 = 0. (15) 
r rmr 8 32 

Considering the initial condition l|13[) . we obtain asymptotical solution of the above equation as follows 

/(r-0) = R + ^r^r-^ + 0(r% 

f'(r^O) = JL r lnr+^r + 0(r 2 ). (16) 



4 



In addition, the minimal surface also implies that /(r) must vanish at a point r — i'q with < r < To < 1, but 
/'(r) is divergent at this point. So that the numerical evaluation will not be valid near this point, and we have to 
find asymptotical solution analytically for r — > ro, 

f(r -> r ) = c(r - r) 1 ' 2 + O((r - rf' 2 ), (17) 

with constant c 2 = 8I(r )G(r )/I'(r ). 

Another aspect which need to be dealt with manually is renormalization of divergence. It is caused by integration in 
Eq. (|10|) when r — > 0. This divergence should be subtracted from numerical result on proper area of string world-sheet. 
Explicitly, using Eqs. I|14|) (|16[) and taking a cut-off A — ► 0, we obtain the divergence as follows 

f R 9 

/ dr ■ r 1/2 fJf' 2 + G(r) = --j— + o A + 0(A 2 ) + other terms independent of A. (18) 

J 2 i / J A 2 L '- > ■ S2R 

Here we include sub-sub leading term which is proportional to X/R such that it will strictly vanish at limit A — ► 0. 
In numerical evaluation, however, the cut-off A is no longer very small. Then this term will play a role for small R. 

According to the above discussions, the Euler-Lagrangian equation (|12J) can be solved numerically. The consistence 
of numerical calculation will be checked via match initial conditions (|13|) and (|17|l . Consequently, the proper area l|l()|) 
of string world can be obtained after subtracting A-dependent terms in Eq. (|18fl . Since difference equation associating 
equation (|12fl is ill-defined for s — > sq, some tricks on numerical evaluation are needed. We put all details on numerical 
evaluation in appendix of this paper. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



We first show some numerical results: 

• Fig. 1 includes some curves of /(r) vs. r, which represents some solutions of Euler-Lagrangian equation i|12|) 
with different initial values R. 
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FIG. 1: Some solutions of Euler-Lagrangian equa- 
tion Jl'it with different initial values. The dash line 
denotes R = 0.32, the dot line denotes R = 2.09 and 
solid line denotes R — 5.06. 
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FIG. 2: Curve of ro vs. R, where ro is zero-point of 
f(r) and R is initial value. 



• Fig. 2 denotes curve of ro vs. R, where ro is zero-point of /(r) and govern integral region in Eq. l|lUfl . i-e., 
< r < r . 

• Fig. 3 and fig. 4 show that the resulted numerical data and related fit curves on proper area A of string world- 
sheet vs. radius R of Wilson loop at large distance and short distance respectively. The unit of area A is taken 
over 3g s N/2 2 / 3 . 

The main conclusions of this paper are included in fig. 3 and fig. 4. In order to precisely describe behavior of 
A ~ In < W(C) > with variation of radius R of Wilson loop, we fit numerical data with several different functions 
and pick up the functions with the smallest x 2 to approach the behaviors of Wilson loop. 
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FIG. 3: Curves of proper area A (in unit 3g a N/2 2/3 ) of 
string world-sheet vs. radius R of Wilson loop at larger 
distance. The diamonds denote numerical results, and 
solid line is the fit function, A = -1.5126 + 0.2941i? + 
1.047? 2 - 2.162 \nR with x 2 = 0.001277. 



FIG. 4: Curves of proper area A (in unit 3g s N/2 2/3 ) 
of string world-sheet vs. radius R of Wilson loop 
at shorter distance. The diamonds denote numeri- 
cal results, and solid line is the fit function, A = 
-1.511 -0.0.829-R + 0.0753 InR with x 2 =3x 10 -5 . 



The resulted functions for long distance (larger R) are listed in Eq. (|19-a|) - (|19-c|l . which fit the numerical data in 
fig. 3. 



A = 


-1.5126 + 


0.2941i? 


f 1.04R 2 - 


2.162 lni? 


x 2 


= 0.0013 


(19-a) 


A = 


-0.7835- 


1.0274R 


f 1.1367E 2 




x 2 


= 0.0046 


(19-b) 


A = 


-0.0464- 


1.717-R + 


1.342i? 2 - 


0.0147? 3 


x 2 


= 0.0017 


(19-c) 



We can see that Eq. (jl9-a[) is the best approach on Wilson loop at large distance. Therefore, it is obvious that 
Wilson loop exhibits area law for large R. In Eq. I|19-af> we also meet linear term and logarithmic term of R which 
are unexpected. These extra terms may be induced by the following reasons: 

1. Due to error bar of numerical evaluation, the subtraction on divergence in Eq. Ijl8(l is not precise. 

2. Other error bar of numerical evaluation. 

3. Due to lack of powerful computers, in our numerical evaluation R is not taken very large. Then resulted Wilson 
loop does not lie in pure confinement phase. These extra may implies some "mixing-phase" . 

However, it is unambiguous that Wilson loop will be dominant by i? 2 -term for very large R. Consequently JV = 1 
SYM is confined at long distance. 

Now let us pay attention to short distance (smaller R). The numerical data and fit curve are shown in fig. 4, and 
fit functions are listed as follows, 

A = -1.511 + 0.829-R+ 0.0753 lni? X 2 = 6.4 x 10~ 4 (20-a) 

A = -1.682 + 1.128i? - 0.131i? 2 X 2 = 0.001 (20-b) 

It is unambiguous that Eq. (|20-a|) is the best fit. Therefore, Wilson loop of M — 1 SYM exhibits logarithmic behavior 
at short distance. It indicates that Af = 1 SYM is asymptotical freedom or approach to Coulomb phase at short 
distance. In addition, we see that the constants in Eq. (|19-a(l and in Eq. I|20-a|) are almost same. Then this constant 
is harmless. It should be associated to renormalization or normalization of Wilson loop. This also indicates fitting in 
Eq. H19-a|) and in Eq. I|20-a|l arc consistent. 

Another interesting issue is that a phase transition occurs when R varies from small ones to large ones. This phase 
transition occurs between the phase of asymptotical freedom and the phase of confinement of non-abelian gauge 
theory. The transition point lies in region 0.5 < R < 2 in unit 3g s N/2 2 ^ 3 . Consequently from Eq. JHJ we can define 
Asym = e 2 / 3 /a' as a scale associating to confinement. Using the results listed in ref . @, Hj, 03 j we have: 
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The masses of glueball and Kaluza-Klein (KK) states scale as 

£ 2/3 

"Vueball ~ ™KK — -— — A SY m/ g S N. (21) 

We usually expect m. g i ue baii ~ i«kk > Asym such that it requires smaller g s ~ 9^ m - In other words, the 
glueball is formed at near perturbative region. Using NSVZ (3 functional for TV = 1 SYM, we have 



^ln^/AlYM + cp 



n 2 N 1 o ln(lnM 2 /A| YM +c ) 1 

a YM 1 — Z 1 9 / a 2 i 



(22) 



where Co is a small constant regularizing divergence at /i = Asym an d fJ- is the scale that glueball or KK states 
are produced (roughly we can set fj, ~ m g i uc baii)- Since the singularity in KST solution is removed through the 
blowing-up of the S 3 of T 1 ' , more precisely, we can include relevant coefficient to Eq. (|21ll . 



3/- 1 / 2 ( So )A S YM/.g s iV, (23) 



where /(so) was defined Eq. (JSJ and s , which denotes the lowest energy detected by certain Wilson loop, is 
zero-point of f(s). Considering radius/energy-scale relation for KST background [l^ . in our nations s ~ /Vmo, 
where the typical scale — Asym consequently it corresponds to Ro ~ 1 and s ~ 0.8. Then defining 
x = tokk/Asym — m/Asym, we achieve the following equation 



3 

6tt- 

l-2 l - 



x^I-y\s /x) ^. lnx2 . (24) 



In x 2 

The solution of the above equation is x ~ 2.6 such that we obtain 

m-glucball — ™kk — 2.6A S ym- (25) 

• The mass of the baryon scales as 

2 /3 

Mb ~ N e — ~ TV A SYM . (26) 
a 

• The gluino condensate of TV = 1 SYM scales as 

< AA > ~ Y^ ~ JVA| YM . (27) 

In QCD wc knew that there is a so-called "confinement scale" , Aqcd — 0.3GeU. Extend to the above results to 
QCD, we have 

"Vuobaii - 2 - 6A qcd ^ 0.8GeU, 

Mbaryon Y C A QCD 0.9GeU, (28) 
< W > ~ A^ CD ~ (0.3GeU) 3 , 

where we use gluino condensate to mimic quark condensate in QCD, but we ignore N c in Eq. (|27|l since those chiral 
quarks are fundamental representation of gauge group. It is surprised that the above results agree with phenomenology 
results of QCD well. It apparently means that pure Kf = 1 SYM mimics some low energy behaviors of QCD even 
without chiral multiplets. 



IV. SUMMARY AND COMMENTS 



Circular Wilson loop of AT — 1 SYM from supergravity KST solution. The Wilson loop exhibits area law at long 
distance and logarithmic law at short distance. Therefore, Af — 1 SYM lies in confinement phase at low energy, and is 
asymptotical freedom or approach to Coulomb phase at high energy. The phase transition occurs when energy varies. 



7 



We also discuss glueball mass, baryon mass and gluino condensate. It is shown that pure Af = 1 SYM mimics some 
low energy behaviors of QCD even without chiral multiplets. 

We would like to make a comment on MN background. In principle, the extension the calculation of this paper 
to MN background is directly. The Euler-Lagrangian equation in MN solution, however, does not possesses initial 
condition like Eq. (|13l) . Instead, it will be / ~ fixed, /' — ► oo at boundary. Since another initial condition at zero-point 
r = tq of f(r) is / — > fixed, /' — * — oo, the numerical evaluation on solution of Euler-Lagrangian is very difficult 
according to criterion in Appendix. The study on this issue, however, is valuable because it can help us to understand 
which phase of J\f = 1 SYM in UV is dual to MN background. Some more discussions on MN background and its 
noncommutative extension are in ref. 201. 



APPENDIX: NUMERICAL EVALUATIONS 



The main difficulty on numerical evaluations are from near-zero point behavior of /(r), i.e., f (r — + ro) — ► — oo. It 
makes Euler-Lagrangian equation l|12|) be stiff for r — > r . If we use Eq. (|13[1 . or more precisely, Eq. (| 1 fc>|) with very small 
r as initial conditions of numerical evaluation. Due to inevitable error on initial condition in numerical evaluation, 
this error will be enlarged rapidly for r — > ro. Consequently we can not obtain a convergent solution. Fortunately, 
since we can analytically obtain solution of Eq. (|12fl for r — * r , we can use Eq. Ill7fl as initial condition of numerical 
evaluation. In this case, the error bar induced by numerical approximation on initial condition is controllable, and it 
will decrease when r is evaluated from ro to zero. So that we obtain the convergent result. 

In our evaluation, we adopt the forth order Runge-Kutta method. We take initial condition at ro — r = 10 -6 and 
interval 2 x 10~ 9 . According to Eq. I|17|l. the error bar on initial condition is about 10 -6 . When we let initial value 
vary 0.1%, the resulted /(0) = R and area only vary 0.0001%. It indicates that the numerical evaluation is indeed 
convergent, and subleading order on initial condition <|17[) can be ignored consistently. 

Another detail is to fix the value of cut-off A defined in Eq. I|18[) . Theoretically, function /'(r) is single- valued in 
region < r < r , and increases from — oo to zero when r varies from r to zero. When string world-sheet lies across 
branes, i.e., in region r < 0, /' (r) is again single-valued and increases from — oo to zero when r varies from — ro to 
zero. It means that f'(r = 0) is maximal value of f'(r) if we extend space to — ro < r < ro. Numerically, however, we 
can not expect to obtain the maximal value of f'(r) at r = exactly and maximal value of f'(r) is not precise zero. 
Then we can take maximal- value point /'(r) as a natural cut-off. In our evaluation, the value of the cut-off A varies 
from 0.0015 to 0.0033 when R = /(0) varies from 0.1 to 5.5. It implies that we indeed consider the second term in 
Eq (|18fl for smaller R in order to achieve high precision fit. 
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